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Abstract 
(N 

1^ ' We compare the cumulative distribution of internal velocity dispersions of galaxy 

^\ . clusters, A^(> cJt,), for a large observational sample to those obtained from a set of N- 

fH \ body simulations that we run for seven COBE^normalized cosmological scenarios. 

Dh; They are: the standard CDM (SCDM) and a tilted (n = 0.85) CDM (TCDM) 

Q '. model, a cold+hot DM (CHDM) model with 9.^ = 0.25, two low-density flat CDM 

(ACDM) models with 9q = 0.3 and 0.5, two open CDM (OCDM) models with 

^ ' f^o = 0.4 and 0.6. The Hubble constant is chosen so that to — 13 Gyrs in all the 

models, while Q.h = 0.02 h? is assumed for the baryon fraction. Clusters identified in 
the simulations are observed in projection so as to reproduce the main observational 
rS ' biases of the real data set. Clusters in the simulations are analysed by applying 

C^ ■ the same algorithm for interlopers removal and velocity dispersion estimate as for 

the reference observational sample. We find that a^ for individual model clusters 
can be largely affected by observational biases, especially for ay^ 600kms~^. The 
resulting effect of A^(> Gy) is rather model dependent: models in which clusters had 
less time to virialize show larger discrepancies between intrinsic (3D) and projected 
distribution of velocity dispersions. From the comparison with real clusters we find 
that both SCDM and TCDM largely overproduce clusters. We verified for TCDM 
that agreement with the observational A^(> ay) requires cts — 0.5. As for the CHDM 
model, it marginally overproduces clusters and requires a somewhat larger a^ value 
than a purely CDM model in order to produce the same cluster abundance. The 
ACDM model with JIq = 0.3 agrees with data, while the open model with f^o = 0.4 
and 0.6 underproduces and marginally overproduces clusters, respectively. 

Key words: galaxies: clusters; cosmology: dark matter, large-scale structure of the 
universe 
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1 Introduction 



The abundance of clusters of galaxies has been recognized in the last few years 
as a crucial constraint for cosmological scenarios of large-scale structure for- 
mation (e.g., Bahcall & Cen 1993; White, Efstathiou & Frenk 1993a); typical 
cluster masses, ~ 5 x lO-'^^^o^"^^© {Hq = 100 /ikms^-'^Mpc^^ is the Hub- 
ble constant), corresponds to a scale length of about lO/i^^Mpc, so that the 
number of clusters with such a mass determines the amplitude of the fluctu- 
ation power spectrum on that scale, once the average matter density is flxed. 
Rather simple analytical arguments, based on the Press & Schechter (1974, 
PS hereafter) approach to the cosmological mass function, demonstrates that 
the cluster abundance is exponentially sensitive to the power spectrum am- 
plitude. This allowed several authors (e.g., Viana & Liddle 1996; Eke, Cole & 
Frenk 1996; Pen 1996) to derive for CDM models tight relationships between 
the r.m.s. fluctuation within a top-hat sphere of 8h~^Mpc radius, ag, and 
the density parameter Qq. Although differing in the details of the derivation, 
they converge to indicate that agfig ' ~ (0.5 — 0.6). Therefore, once the 

value of Qq (and of the cosmological constant term) is chosen, requiring at the 
same time to satisfy the as constraints, coming at small scales from the cluster 
abundance and at large scales from measurements of CMB anisotropics (see, 
e.g.. Bond & Jaffe 1996; Lineweaver et al. 1996), represents a stringent test 
for the shape of the power spectrum. The most famous example is probably 
represented by the standard CDM model, that, once normalized to match the 
large-scale CMB anisotropics, overproduces clusters by at least one order of 
magnitude. 

It is however clear that, in order to fully exploit the potential of the cluster 
mass function as a cosmological test, one has to be sure about (a) the reliability 
of the PS approach, and (b) the accuracy of cluster mass measurements. As 
for the PS approach, its accuracy has been verifled against a variety of N- 
body simulations (e.g.. White et al. 1993a; Lacey & Cole 1994; Borgani et al. 
1997a). As for cluster masses, methods for their determination are based on 
virial analysis of internal galaxy velocities. X-ray temperature observations 
and gravitational lensing of background galaxies. At present, the above three 
methods lead often to discordant results for the same cluster (e.g., David, 
Jones & Forman 1995, Miralda-Escude & Babul 1995; Wu & Fang 1997, and 
references therein). 

Masses based on X-ray data are usually considered as the most reliable (e.g., 
Evrard, Metzler & Navarro 1996, and references therein) and are widely em- 
ployed to constrain cosmological models (see, e.g., Oukbir & Blanchard 1996, 



and references therein). However, passing from temperature measurements to 
cluster masses relies on assumptions, like spherical symmetry and hydrostatic 
equilibrium, which may not be met in realistic cases. Observations of substruc- 
tures in the temperature patterns (e.g., Mohr et al. 1995; Henry & Briel 1995; 
Buote & Tsai 1996) and in the internal galaxy distribution (e.g.. Dressier & 
Schechtman 1988; Escalera et al. 1994; Bird 1995; Crone, Evrard & Richstone 
1996) indicates that clusters may not be such completely relaxed structures, 
thus casting doubts on the robustness of mass determinations based on hydro- 
static equilibrium (cf. Balland & Blanchard 1996). Furthermore, the possible 
presence of pressure supports of non-thermal origin, like intra-cluster mag- 
netic fields (e.g., Loeb & Mao 1994; Ensslin et al. 1996) could also lead to a 
mass underestimate. 

As for gravitational lensing, there are at present only a few clusters at moder- 
ate or high redshift whose mass is estimated from this method. Furthermore, 
weak lensing observations are more reliable in providing the shape of the in- 
ternal mass distribution, rather than mass estimate (e.g. Squires & Kaiser 
1996; Seitz & Schneider 1996). Measurements of the overall cluster masses are 
based on strong lensing, whose application is however limited to very central 
regions. In any case, all these estimates are rather dependent on the details of 
the lens model (e.g. Bartelmann 1995). 

Optical virial mass estimates have also been attempted by several authors 
(e.g., Biviano et al. 1993). However, even in this case, passing from the internal 
velocity dispersion of galaxies to the cluster mass requires suitable assump- 
tions about the degree of virialization of the cluster, the nature of the galaxy 
orbits, the relation between mass and galaxy distribution profiles (e.g., Mer- 
ritt 1987) and, possibly, the presence of substructures (Bird 1995). A detailed 
analysis of biases affecting cluster mass estimates from galaxy velocity disper- 
sions have been recently pursued by Frenk et al. (1996). From the analysis of 
numerical simulations they concluded that virial methods can underestimate 
cluster masses by up to a factor five. 

Due to such uncertainties, a preferable procedure is to resort directly on the 
line-of-sight velocity dispersion (o"„, hereafter), which is the observable quan- 
tity, instead of on the cluster mass, as a diagnostic for cosmological models. 
Cluster velocity dispersions are rather easily obtainable from numerical N- 
body simulations and, in fact, their distribution has been used as a constraint 
for dark matter scenarios (Frenk et al. 1990; Bartlett & Silk 1993; Jing & Fang 
1994; Crone & Geller 1994). 

However, several observational biases may well affect a^ estimates, for instance 
connected to projection effects (e.g.. White et al. 1990; van Haarlem, Frenk 
& White 1997) and to the limited number of galaxies usually available to 
trace the internal cluster dynamics (see, e.g., Mazure et al. 1996 and Fadda 



et al. 1996, for recent detailed discussions about robust methods to estimate 
a^). Moreover, cluster member assignment in projection, anisotropy of galaxy 
orbits, cluster asphericity, infalling galaxies and presence of substructures can 
significantly affect any measurement of a^ and must be taken into account to 
reliably compare data and simulations. 

A further problem in the determination of the a^ distribution concerns the 
availability of an observational sample which has to satisfy well defined com- 
pleteness criteria, to be implemented also in the simulated sample. All the 
up-to-date available estimates of the cumulative velocity dispersion function 
(CVDF hereafter), A^(> a^), are based on samples that are complete with 
respect to cluster richness (Girardi et al. 1993; Zabludoff et al. 1993; Collins 
et al. 1995; Mazure et al. 1996, M96 hereafter). However, since the relation 
between richness and a^ is intrinsically rather broad, all samples that are com- 
plete in richness are not complete in a^ (see Figure 8 of M96, see also Girardi 
et al. 1993). Instead, they are systematically biased toward low a^^s. For in- 
stance, the ESO Nearby Cluster Survey (ENACS; Katgert et al. 1996), which 
includes all the clusters with i? > 1, is complete in velocity dispersion only for 
cr„^ 800 km s~^. In spite of all these problems, there is a fair agreement among 
observational distributions coming from different studies, at least within the 
completeness limits of the corresponding samples (see, e.g., Fadda et al. 1996; 
F96 hereafter). It is however clear that such observational limitations must be 
taken into account in order to reliably compare real data to simulations. 

This paper is devoted to a close comparison of the cr^ distribution from the 
sample analysed by F96 and from an extended set of particle-particle-particle- 
mesh (P3M) simulations based on seven different cosmological models. The 
original code on which the simulations are based is the adaptive P3M one 
provided by Couchman (1991). 

As for the F96 sample, on which our analysis is based, it includes 153 Abell- 
ACO (Abell, Corwin & Olowin 1989) clusters with redshift z < 0.15. Data 
were taken from the published literature and from the ESO Nearby Abell 
Clusters Survey (ENACS; Katgert et al. 1996). All clusters were selected so 
as to have at least 30 galaxies with measured redshift within the cluster field 
and a significant peak in redshift space (> 99% c.l.; see F96 for a detailed 
description of the sample). The inclusion of clusters of lower richness with 
respect to ENACS pushes downwards the completeness limit for a^. As a 
result, F96 estimate their CVDF to be complete above (Ty^um — 650kms~^. 

Clusters extracted from the simulations are treated so as to reproduce the main 
features of the observational analysis. Firstly, they are observed in projection 
with the same aperture radii and number of galaxies with measured redshift 
as in the F96 sample. Then, we applied an algorithm for removing interlopers 
and estimating cr„ which is the same as that of F96. 



The aim of our analysis is twofold. From the one hand, the availability in the 
simulations of the whole 3D information about the cluster internal dynamics 
provides us with a test for the robustness of the N[> ay) estimate. From the 
other hand, the possibility of analysing an extended set of simulations based, 
on different initial spectra of fluctuations, allows us to put useful constraints 
on the parameter space of cosmological models. 

The plan of the paper is as follows. 

We describe in Section 2 the simulations and the models that we considered. 
In Section 3 we outline the procedure for cluster identification in 3D and 
how we simulate cluster observations to reproduce the features of the F96 
sample. After re-analysing the F96 sample, in Section 4 we show in detail how 
the algorithm, which estimates a^ from projected data, works and discuss 
the effects of observational biases on the A^(> o"^) for the different models. 
Afterwards, we compare the CVDF for data and simulations and discuss the 
resulting constraints on cosmological scenarios. A summary of the main results 
and our conclusions are deserved to Section 5. 



2 Models and simulations 



As already mentioned the standard version of the CDM model (SCDM) turns 
out to largely overproduce the abundance of galaxy clusters, once it is normal- 
ized on large scales to match the measured CMB anisotropics. Therefore, we 
decided to simulate SCDM, as a kind of reference model, as well as six more 
models, which correspond to different ways of improving SCDM. Such models 
can be divided into three main categories as follows. 

(a) Two r^o = 1 models, namely a purely CDM model having a "tilted" 
primordial spectral index (TCDM), n = 0.85, and a Cold+Hot Dark Matter 
(CHDM) model with Q^ = 0.25 for the density parameter contributed by 
one species of massive neutrinos with m,y ~ 6.25eV. 

(b) Two flat low-density CDM models, characterized by {Qq, h) = (0.3, 0.73) 
(ACDM0.3) and ino,h) = (0.5,0.63) (ACDM0.5), respectively, with the flat- 
ness provided by the cosmological constant term Q/^ = 1 — fig. For both 
models, the Hubble parameter h is tuned in such a way to give tg — 13 
Gyrs, so as to be consistent with the f^o = 1 models. 

(c) Two open CDM models, with (fio, h) = (0.4, 0.59) (OCDM0.4), and (l^o, h) -- 
(0.6,0.55) (OCDMo.e), with the same criterion as before for the choice of 
the Hubble parameter. 

A further way to suppress cluster abundance in CDM models, that we do 
not explore here, would be to abandon scenario of random phase adiabatic 



fluctuations. For instance, a (^o = 1 CDM cosmology with topological defects 
can easily provide o-g — 0.5 (e.g., Pen, Spergel & Turok 1994). 

For the CDM models we adopted the transfer function by Bardeen et al. 
(1986), with shape parameter F = ilohexp{—ilb — (h/O.bY^'^ilb/^o), corrected 
to account for a non negligible baryon contribution (Sugiyama 1995). The 
baryon fraction has been chosen to be Qb = 0.02h~^ for all the models. This 
value corresponds to the 95% upper limit from primordial nucleosynthesis 
predictions (see, e.g., Copi, Schramm & Turner 1995; see Buries & Tytler 
1996, for higher Q^ from observations of low deuterium abundance in high 
redshift systems; see, however. Fields et al. 1996 for a lower Qb prediction). As 
for the CHDM power spectrum, it has been explicitely computed by following 
the linear evolution of the matter and radiation fluids through the epoch of 
matter-radiation equality and recombination epochs, down to the redshifts 
relevant for large-scale structure formation. We normalized all the spectra 
to match the four year COBE data (e.g., Bennett et al. 1996), following the 
recipe provided by White & Scott (1996). The resulting values for the r.m.s. 
fluctuation amplitude within a top-hat sphere of 8 h~^Mpc radius, ag, are 
reported in Table 1, where also the model parameters are specifled. 

Also reported in the last column of Table 1 are the erg values predicted by 
the Eke et al. (1996) [E96 hereafter; see their eqs. (6.1) and (6.2)] to re- 
produce the distribution of X-ray cluster temperatures by Henry & Arnaud 
(1991), based on the Press & Schechter (1974) approach and the assumptions 
of isothermal gas distribution and hydrostatic equilibrium. Other approaches 
have been derived by different authors, who provided slightly different results. 
For instance, Viana & Liddle (1996) found that ag — 0.6 is required for Oq = 1 
when normalizing the temperature function only at Tx = 7 keV. Pen (1996) 
pointed out that the E96 scaling provides too small erg values ii Qq < 1. For 
Q\ = 0.65 he found the ag value by E96 to be underestimated by about 17%. 
In general, such ocnstraints from X-ray data agree with the cluster abun- 
dance as inferred, although with larger uncertainties, from the frequency of 
large-separation lenses (see, e.g., Kochanek 1995). Therefore, we regard the 
numbers reported in the last column of Table 1 more as guidlines than as 
stringent constraints. Having this in mind, we note that our models are not in 
general tuned so as to reproduce such predictions. Indeed, our purpose here 
is to verify whether a careful analysis of cluster velocity dispersions leads to 
the same conclusions as the above analytical approaches, rather than picking 
up the best-fltting cosmological scenario. 

The reported o-g value for TCDM assumes a vanishing tensor (gravitational 
wave) contribution to the CMB anisotropy. In addition to this, we also con- 
sider two more outputs at cg = 0.67 and ag = 0.51. The flrst value corresponds 
to assuming T/S = 7(1 — n) for the ratio between tensor and scalar contribu- 
tions to the CMB anisotropy, as predictied by power-law inflation (see, e.g.. 



Table 1 

The model parameters. Column 2: the density parameter 0,q; Column 3: the cos- 

mological constant term Q\; Column 4: the Hubble parameter h; Column 5: the 4 

year COBE predicted linear r.m.s. fluctuation amplitude at 8/i^^Mpc crgj Column 

6: the initial redshift; Column 7: The cxg value predicted by Eke et al. (1996; see 

text). 



Model 


Oo 


^A 


h 


(^8,COBE 


Zi 


(^8,E96 


SCDM 


1.0 


0.0 


0.50 


1.20 


35 


0.52 ± 0.04 


TCDM 


1.0 


0.0 


0.50 


0.85 


25 


0.52 ± 0.04 


CHDM 


1.0 


0.0 


0.50 


0.77 


20 


0.52 ± 0.04 


ACDMo.3 


0.3 


0.7 


0.73 


1.10 


40 


0.94 ± 0.07 


ACDMn 5 


0.5 


0.5 


0.63 


1.22 


40 


0.73 ± 0.06 



OCDMo.4 0.4 0.0 0.59 0.58 35 0.77 ±0.06 
OCDMo.e 0.6 0.0 0.55 0.90 35 0.64 ± 0.05 



Crittenden et al. 1993). The second value is consistent with the observational 
constraint reported in column 7. 

We simulate the development of non-linear gravitational clustering by resort- 
ing to the adaptive P3M code developed by Couchman (1991). For each model 
we run three realizations within a box of size L = 150 /i~^Mpc, while only one 
realization is run for CHDM. The evolution of the density field is traced by 
following the trajectories of 128^ cold particles. In additions to these, for the 
CHDM model we also put 2 x 128^ hot particles, in order to sample the neutrino 
phase-space. Therefore, for all the purely CDM models the particle mass is 
4.5 X lO^^Qoh^^MQ, while cold and hot particle masses for the CHDM simula- 
tion are about 3.4 x IO^^H^^Mq and 5.6 x IO^^H^^Mq, respectively. Therefore, 
a typical cluster of mass ~ IO^'^H'^Mq is resolved with more that 200 r2(^^ 
particles, thus ensuring an adequate mass resolution. 

Initial conditions are realized by generating a random realization of the lin- 
ear gravitational potential on a 128^ grid. Particles are moved from their 
initial grid position q according to the Zel'dovich approximation x = q — 
J2]^ k0fce~*"^'^ with comoving peculiar velocity u = x = —a J2]^ cifck^j^e^*"'^'^, 
where a^ = (ilog(5]^/(iloga. Accordingly, a^ = const for purely CDM models 



(ttfc = 1 if fio = 1), while it has to be exphcitely computed for the CHDM 
model in order to account for the effect of residual neutrino free streaming 
on the fluctuation growth (see also Ma 1996; Klypin, Nolthenius & Primack 
1996). 

As for the hot particles in the CHDM simulation, in addition to the gravita- 
tional velocities, they were given also a thermal velocity which is randomly 
taken from the Fermi-Dirac distribution 



fpoiv) oc 



v^ 



exp [m^v/kBT^izi)] + 1 



Here rriu ~ 6.25eV is the neutrino mass corresponding to Q^, = 0.25, ks is the 
Boltzman constant and T„{zi) = (1 + Zi)T^fi, being T^^ = {A/liy^^TcMB,o — 
1.95 K the present-day temperature of the relic neutrino background. Each 
pair of hot particles, which are initially located on the same grid position, are 
assigned thermal velocities having the same modulus but opposite directions, 
so as to ensure local momentum conservation. 

As for the long range force, it is computed on a 128^ grid, while the short- 
range force is softened at scales smaller than e = 0.1(L/128). Simulations are 
started from an initial redshift Zi corresponding to the epoch at which a = 0.2 
for the linearly estimated r.m.s. fluctuation on the grid. The resulting values 
of Zi are listed in column 6 of Table 1. The integration variable has been taken 
to be p = a^/^, being a = {I + z)/{l + z,i) the expansion factor. We adopted a 
constant step size Ap ^ 0.35 for all the models, except for CHDM, for which 
we take Ap ~ 0.3 to allow for a slightly more accurate integration. Thanks 
to their dynamical and mass resolution, our simulations are well suited to 
follow the internal dynamics of galaxy clusters on the scales relevant for the 
ay estimate. 



3 Construction of the cluster samples 

3. 1 Cluster identification in 3D 



Our method to identify clusters in the simulation box is based on the friend- 
of-friend (FOF) algorithm. The list of candidate clusters is constructed by 
finding groups of CDM particles using a linking length h = 0.2 (in units of 
the mean particle separation; see, e.g., Frenk et al. 1990; Lacey & Cole 1994). 
In order to approach a more observational selection procedure, for each FOF 
group we estimate the center of mass and draw around it a sphere having 
radius equal to the Abell one, tm = 1.5/i^^Mpc. The center of mass of all 



the cold particles falling within this sphere is then computed and used as the 
starting point for the next iteration. The sphere is moved around until we get 
convergence for its mass and position. We always find that few (~ 5) iterations 
are required to obtain a stable result. Cold particles are used to trace their 
internal velocity dispersion. When the distance between the centers of two 
clusters is smaller than 2rAb, the less massive cluster is removed from the list. 
We find that this situation occurs in most cases when a small group is located 
at the outskirts of a massive object. 

For the above choice of the FOF linking length, the group-finding algorithm 
picks up structures with p/p ~ 180 Qq for the average internal overdensity, 
thus very close to the virialization overdensity for f2o = 1. If fig < 1, smaller 
b values are required to pick up virialized structures (see, e.g., E96); it turns 
out that b ~ 0.15 is required in the most extreme case of the ACDMq.s model. 
However, since FOF groups are only used as initial points from where to start 
looking for clusters, one expects the final cluster list not to be very sensitive 
to the choice of b. We verified that final results are insensitive to variations of 
b in the range 0.15-0.2. 



3.2 Observing simulated clusters 



The list of clusters identified in 3D is used as the starting point to obtain a 
sample which resembles as close as possible the features of realistic cluster 
observations. The procedure to pass from the 3D to the projected estimate of 
the cluster velocity dispersion is realized according to the following steps. 

(1) Each 3D cluster is observed by selecting all the particles contained within 
a cylinder of fixed aperture radius, r^ = 1.5/i^^Mpc, which extends for 
±8, 000 km s~^ in line-of-sight velocity from the cluster center, by allowing 
for periodi boundary conditions. Each cluster is observed three times along 
the directions of the three coordinate axes. Galaxy velocity distributions for 
observed clusters do not extend in general outside ±4000 km s^^ from the 
mean cluster velocity (see, e.g., Zabludoff, Huchra & Geller 1990). Hence, 
the extension of ±8000 km s^^ in l.o.s of our selected cylinders is sufficiently 
large to allow the density peak finding algorithm to retrieve cluster peaks. 
In particular, we avoid spurious effects at the border of the velocity range, 
that may occur when a peak of a close cluster is sharply truncated there. We 
verified that this happens few times if a smaller extension of ±4000 km s~^ 
is considered. 

(2) In the observational case the sampling density of cluster galaxies is by 
far much smaller than that allowed in the simulations and improves with 
cluster richness (as we verified by using the data sample analysed by F96). 
Furthermore, an aperture radius as large as 1.5 h^^Mpc is in general ap- 



propriate only for very rich clusters (see, e.g., Katgert et al. 1996), poor 
clusters and groups requiring r^ ~ 0.5 h~^Mpc (e.g., Dell' Antonio, Geller & 
Fabricant 1995). Therefore, observing the latter at a much larger aperture 
one may include a nearby close cluster, whose presence would pollute the 
velocity dispersion estimate. To account for this, we group clusters of the 
observational sample (F96) into four equally spaced bins in o"„, from to 
1300 km s~^. Then, we assign to a model cluster with velocity dispersion a^ 
the aperture radius r^ and the number of galaxies Ng of a real cluster which 
is randomly selected between those belonging to the same a^ bin. In this 
way, we generate a r^-cr^ and a Ng-a^ correlation, having the same shape 
and dispersion as for the observational sample. 
(3) Galaxy membership to clusters is assigned by following the same proce- 
dure of F96. 

Starting from the discrete distribution of galaxies in velocity space within 
each cylinder, a continuous distribution is obtained by applying the adaptive 
kernel method (see Pisani 1993, 1996 and references therein). This method is 
based on convolving the discrete distribution with a Gaussian kernel, whose 
r.m.s. amplitude is chosen point by point so as to minimize the difference 
between the original distribution and its reconstructed continuous represen- 
tation. Peaks of this continuous distribution which are significant at > 99% 
level are then identified with clusters. In the case of two or more peaks, we 
retain only the one closer to the 3D cluster position on which the cylinder 
is centered, so as to avoid multiple counting of the same object. 

Then we applied the "shifting gapper" method which is an iterative pro- 
cedure based on combining projected position and velocity information for 
each galaxy. A galaxy is considered as interloper if it is separated by more 
than 1000 km s~^ from the central body of the velocity distribution of galax- 
ies lying at the same distance (within a bin of 0.4 h~^Mpc width) from the 
cluster center. Therefore, this method accounts for the possibility that the 
velocity dispersion strongly depends on the clustercentric distance, as seen 
in several cases for observed clusters (see, e.g., F96 and Girardi et al. 1996, 
G96 hereafter). 

Once the cluster membership has been assigned, the corresponding a^ 
is estimated by using the robust estimator described by Beers, Flynn & 
Gebhardt (1990). 

As in F96, we discarded from the sample those clusters whose a^ has a 
bootstrap error larger than 150kms~^. A visual inspection of the velocity 
dispersion profile of such clusters systematically reveals that these large 
errors occur when o"„ widely oscillates at large radii, instead of smoothly 
converging to a constant value. This could be the signature for an insufficient 
sampling or for the presence of a heavy contamination due to unremoved 
interlopers, whose effect is that of boosting o"„ for an otherwise low velocity 
dispersion group. 
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Although this procedure is designed so as to reproduce the observational situa- 
tion as close as possible, it is clear that several aspects have not been included, 
due to intrinsic limitations of our simulations. In the following we discuss these 
limitations and comment their effects on the final results. 

(a) The projected richness of optically selected clusters can be in general a bad 
predictor of the real 3D richness (see, e.g., Frenk et al. 1990). For instance, 
van Haarlem et al. (1997) pointed out that (i) about one third of i? > 1 
Abell clusters identified in Dq = 1 CDM N-body simulations can arise from 
superposition of intrinsically poor clusters, and (ii) 30% of intrinsically rich 
clusters are missed and classified as poor clusters because of fluctuations 
of the background galaxy counts. On the contrary, Cen (1996) addressed 
the same issue by analysing an Qq = 0.4 open CDM simulation and argued 
against a substantial richness contamination by projection. Whether such 
a disagreement is due to the different nature of the simulated models (in 
an open Universe clusters are smaller, more concentrated structures than 
for r^o = 1) or to differences in the simulations, we regard this as an open 
question. As far as our comparison with observational data is concerned, 
we point out that misclassified clusters ought not to represent a serious 
problem. In fact, even though a cluster can be misclassified in projection, 
it can be recognized as a superposition of poor clusters when observed in 
redshift space (e.g. in the case of A151 and A367, see F96). As for the missed 
rich clusters, M96 showed that the effect of this bias is that of reducing the 
range of completeness of the observed 0"„ distribution thanks to the existence 
of a broad, but well defined, correlation between the cluster richness in two 
and three dimensions (see Figure 7 in their paper). Since our reference 
CVDF is estimated to be complete down to (Jy^im. — 650 km s~^ (see F96) 
we account for this bias by comparing A^(> cr^) for data and simulations 
only at a^ > a^^um (see § 4.2 below). This avoids any need to define the 
richness for model clusters, which necessarily passes through a problematic 
galaxy identification in the simulations. 

(b) In their analysis, F96 removed the contribution of late-type galaxies, 
which is expected to be more affected by interlopers, when their velocity 
dispersion is significantly larger than that of early-type galaxies. However, 
since any distinction between early- and late-type galaxies is not allowed in 
our simulations, the analysis of the observational sample has been repeated 
by including galaxy of any morphology (see § 4.1). We expect this to turn 
into an increase of the a^ estimates for the combined effects of a larger 
number of interlopers and of the larger velocity dispersion of spirals with 
respect to ellipticals (cf. also Stein 1996). 

Furthermore, F96 also attempted to remove the contribution of velocity 
gradients, possibly produced by cluster substructures, cluster rotation or 
by the interaction with surrounding structures, like filaments and nearby 
clusters. Since such effects are already accounted for in the simulations, we 
have chosen not to correct for them in the reanalysis of the observational 
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sample (see §4.1). 
(c) Our comparison between data and simulations assumes that galaxies are 
fair tracers of the internal cluster dynamics, that is, no velocity bias is in- 
troduced. Whether this represents a reasonable approximation is a matter 
of debate. Carlberg (1994) found in its purely gravitational simulation that 
galaxy halos suffer for a substantial velocity bias. This result has been con- 
firmed by Evrard, Summers & Davis (1994) and by Frenk et al. (1996), who 
also included gas dynamics. On the contrary, Katz, Hernquist & Weinberg 
(1995) and Navarro, Frenk & White (1995) came to a different conclusion, 
claiming for no substantial velocity bias in their cluster simulations. (Note 
that the mass resolution of our simulation does not allow to address this 
issue.) As for observations, comparisons of galaxy velocity dispersions and 
X-ray temperatures in clusters (Edge & Stewart 1991; Lubin & Bahcall 
1993; G96; Lubin et al. 1996) indicate that the former provide a fair tracer 
of the dark matter velocity dispersion (see, however. Bird, Mushotzki & 
Metzler 1995). In the following we will show explicitely only for the SCDM 
simulations the effect of introducing the velocity bias found by Evrard et 
al. (1994) for their cluster simulations based on this same model. 



4 Results 



4-1 The Observational Distribution of Velocity Dispersions 



In order to allow for a more homogenous comparison with simulation results 
(see point (b) of §3.2), we used as a starting point the F96 sample in an 
intermediate step of its analysis, that is before applying removing late-type 
spirals and correcting for velocity gradients (see § 2.2 and § 3 of F96). More- 
over, we considered only galaxies observed within a maximum aperture radius 
of 1.5/i~^Mpc. In the case of multipeaked clusters (see § 4 in F96), we se- 
lected only the most significant peak. We computed a^ for each cluster and 
the resulting CVDF by following the same procedure as F96. In particular, 
we weighted each cluster according to the richness-class distribution of the 
Edinburgh-Durham Cluster Catalogue (EDCC; Lumsden et al. 1992) in order 
to account for the volume-incompleteness of the observational sample (see § 5 
of F96). 

We compare in Figure 1 the original A^(> a^) from F96 (squares) to the one 
we calculated here (circles). Also plotted is the CVDF obtained by M96 for the 
EN ACS clusters (continuous line). As already pointed out in the Introduction, 
the ENACS sample is incomplete for a^^ 800 km s^^. This explains the bend- 
ing of its N[> a^) in this velocity range. As for the F96 sample, it is remark- 
able the agreement with the original A^(> cr„) estimate. This indicates that, at 
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Fig. 1. The cumulative velocity dispersion function for the observational sample as 
originally published by F96 (squares) and re-estimated here (circles). The continu- 
ous line represents the determination by Mazure et al. (1996) for ENACS clusters. 
Errorbars are the la uncertainties obtained by adding in quadrature errors coming 
from the Poissonian statistics and from the bootstrap error in the determination of 
individual cluster velocity dispersions. 

least from a statistical point of view, velocity gradients and galaxy morphol- 
ogy play a marginal role. In turn, these results are very close to those by M96 
in the range of completeness of the ENACS. The plotted errorbars (reported 
only in one case for sake of clarity) represent the la uncertainties obtained by 
summing in quadrature the Poissonian errors and the bootstrap errors in the 
determination of ay for each individual cluster. Note that Poissonian errors 
are due to the fact that we sample the intrinsic A^(> a^) distribution with 
a finite number of clusters, while bootstrap errors are connected to the finite 
sampling of the velocity dispersion within individual clusters. Therefore, such 
errors have independent origin and must both be considered. 
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In the following, our CVDF will be used for the comparison with simulation 
results. 



4.2 Testing the method 



We will test here the reliability of the method for cluster membership assign- 
ment, that we described in the previous section. We will also compare the 
resulting o"„ estimates to the intrinsic ones, as provided by using the whole 
3D information from the simulation particle distribution. After analysing how 
projection and sampling biases affect the determination of the velocity disper- 
sion, we will also show their overall imprint on A^(> a^). 

We plot in Figure 2 the result of introducing observational biases on four 
clusters selected from the TCDM model with cr^ = 0.67. The first three clus- 
ters from the left have progressively larger a^ values, while the fourth column 
shows the case of a cluster which is rejected on the ground of its too large o"„ 
bootstrap error [> 150kms~^; cf. point (3) in the previous section]. Panels in 
the first line show how the clusters appear in projection, with open circles in- 
dicating the rejected interlopers and filled circles for the "galaxies" recognized 
as genuine members. We note that the algorithm for membership assignment 
is rather efficient in identifying interlopers preferentially at the outskirts of the 
clusters, while genuine members are more concentrated to define the cluster 
shape. This is also confirmed by the panels in the second line, which show the 
real-space distribution along the line of sight for accepted (filled histogram) 
and rejected (open histogram) members. In general, true members of the first 
three clusters are correctly recognized to lie very close to the cluster cen- 
ter, with few exceptions (e.g., note the group of seven galaxies, located at 
~ 20/i~^Mpc from the center of the second cluster). As for the fourth clus- 
ter, a non-negligible fraction of galaxies recognized as genuine members are 
instead interlopers, which lie at a distance ^ 20 h~^Mpc from the cluster. 

Panels in the third line show the redshift-space distribution from which o"„ 
is actually calculated, while the fourth line shows the projected integrated 
velocity dispersion profile. In each panel the horizontal dashed line represents 
the intrinsic value of a^, which is estimated from all the simulation particles 
lying within 1.5/i~^Mpc from the cluster center. Errorbars are the la scatter 
from 1000 bootstrap resamplings. As for the first two clusters, it is remarkable 
how well the correct a^ is recovered at the scales where it keeps fiat (see F96 
and G96 for discussions about the fiatness of profiles in the external cluster 
region). As for the third cluster, despite the fact that the a^ profile fiattens 
at scales ^ 1 /i~^Mpc, it does not converge to the correct value. Note that 
the discrepancy of ~ 200kms~^ exists at a high confidence level (> 3a) and, 
therefore, can not be ascribed to sampling uncertainties. This means that 
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Fig. 2. Results of the application of the cluster member assignment and a^ estimate 
on four clusters extracted from the TCDM simulations with erg = 0.67. Panels in 
the first line show clusters in projection. Filled and open circles refer to "galaxies" 
recognized as genuine members and interlopers, respectively. Panels in the second 
and third lines show the histograms for the line~of-sight galaxy distribution toward 
the cluster centers in real and velocity space, respectively. Open and filled histograms 
are for the distributions of interlopers and true cluster members. Panels in the 
fourth line are for the integral velocity dispersion profiles. The horizontal dashed 
line corresponds to the intrinsic ay value, computed at 1.5h~^M.pc. Errorbars are 
for the r.m.s. scatter over 1000 bootstrap resamplings. For each cluster the point at 
the smallest radius represents a^ as computed for the first six innermost galaxies. 

a better recovering of the intrinsic a^ can hardly be attained in this case 
with a denser galaxy sampling. The situation is quite different for the fourth 
cluster. In this case, the discrepancy is as large as ~ 400kms~^, but it is 
not significant, as a consequence of the large errorbars. This example shows 
the effectiveness of eliminating those clusters with large (T„ errors, which in 
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general corresponds to rather small and loose structures, whose a^ can be 
heavely boosted by sampling uncertainties. 

Figures 3-5 are devoted to the statistical description of what can be under- 
stood from measurements of cluster internal velocity dispersion. Results are 
reported for all the nine simulation outputs that have been considered (for 
each model, plotted results refer to only one realization). The panels in the 
first line show the relation between the intrinsic (7v,zD) estimated by using all 
the simulation particles within 1.5/i~^Mpc, and the total cluster mass within 
the same radius (as for CHDM, both quantities refer only to the cold par- 
ticles). We note that a well defined correlation always exists, which reflects 
the condition of virial equilibrium characterizing most of the clusters. Few 
clusters detach significantly from this correlation and all have velocity disper- 
sions larger than the virialization one. Such exceptions are even more rare for 
those models, like the OCDM ones, in which the cluster particles have spent 
more time within the collapsed structure and, therefore, have had more time 
to reach virial equilibrium. 

Note that a scatter of a factor two in mass at a fixed a^ value is not rare, 
especially at low a^ values. However, one would conclude that in general rather 
reliable cluster mass determinations would be possible from measurements of 
the internal velocity dispersion. The situation is rather different in realistic 
cases, when observational limitations are introduced. This is shown in the 
panels of the second line, where we plot the relation between the projected 0"^,, 
estimated according to the procedure described in Section 3, and the intrinsic 
(Tj,. Although most of the clusters show discrepancies of about 100-200 km s~^ 
between observed and true velocity dispersions, there are cases in which the 
difference is as large as ~ 500kms~^. 

The consequence of the large scatter in the (Jv,proj^(^v,?,D on A^(> a^) is plotted 
in the panels of the third line, where we compare the CVDF for the intrin- 
sic (solid lines) and observed (dashed lines) velocity dispersions. Errorbars, 
which are reported only for the "observational" CVDF, represent the Pois- 
sonian uncertainties. As expected, the overall effect is that of increasing the 
CVDF, especially in the large (T„ tail. This result is qualitatively similar to 
that found by van Haarlem et al. (1997). A closer comparison with his anal- 
ysis is however rather difficult, since he considered simulations only of the 
SCDM model with erg = 0.63, by adopting a lower mass resolution, although 
within a larger simulation box. Furthermore, he also used different procedures 
to remove observational biases in the a^ estimate. 

We point out that all the results we reported here are based on clusters whose 
?)D velocity dispersion is larger than aiim,3D = 400kms~^. Since we require 
that simulated cluster samples be complete above aiim,proj = 650kms~^ (in- 
dicated with the horizontal dashed line), we should keep aum^sD small enough 
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Fig. 3. Effect of the observational biases on the a^ estimates for the three TCDM out- 
puts. Panels in the first line show the relation between intrinsic velocity dispersion 
and cluster masses, both estimated at 1.5 /i~^Mpc. Panels in the second line show the 
scatter between intrinsic and observed cluster velocity dispersions. The horizontal 
dashed line is for the a^ completeness limit of the F96 sample, cr^^um = 650 km s^^. 
Panels in the third line are for the intrinsic (solid lines) and observational (dashed 
lines) CVDF. Errorbars represent the la Poissonian uncertainty. 
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Fig. 4. The same as in Figure 3, but for SCDM, ACDMq.s and ACDM0.5 simulations. 



that no clusters below this limit would have their a^^proj increased above the 
observational completeness limit. On the other hand, the number of selected 
clusters rapidly increases as (Jum^^D decreases, so as to hardly keep the amount 
of data to be analysed to a manageable size. We verified that the incomplete- 
ness induced in the simulated samples by taking (Jum,^D = 400 km s^^ is always 
negligible (note that only very few clusters have a^^proj > 650kms~^ at the 
smallest o"^,3d)- Indeed, for the TCDM model with erg = 0.67 we plot the ob- 



served A^(> a^) obtained from cTiim,3D = 400kms~^ (short-dashed curve) and 
Ciim,3D = 200kms~^ (long-dashed curve). The small difference between these 
two cases confirms the reliability of cutting at (Jiim,zD = 400kms~^. 

Note that the difference with respect to the intrinsic A^(> a^) is rather model 
dependent being in general smaller for those models whose clusters have had 
more time to virialize; an excellent recovering of the 3D cumulative distribu- 
tion is indeed found for SCDM and ACDM0.5, which have the largest ag values; 
instead, larger differences exist for TCDM as lower ag values are considered. 
Also note how models whose intrinsic A^(> o-^,) are comparable, like ACDM0.3 
and TCDM with a^ = 0.51 are affected in a rather different way by observa- 
tional biases. The CVDF for ACDM0.3 is almost unaffected at ay^ 900kms~^, 
while it only acquires a high cr^ tail up to cr^ ~ 1100 km s~^. On the other hand, 
the low-CTg TCDM significantly increases its CVDF over the whole a^ range. 
This is due to the fact that ACDM0.3 clusters are small and rather isolated 
structures, which already undergone virialization. As a consequence, the ef- 
fect of interlopers is in general smaller than for TCDM clusters, which are 
expected to be more extended objects (like in any Dq = 1 model) character- 
ized by the presence of substructures and continuous merging of surrounding 
clumps. This is explicitely shown in Figure 6, where we plot the structure and 
velocity dispersion profiles for two clusters extracted from a ACDM0.3 and a 
TCDM (o"g = 0.51) simulation. Left and central panels show in projection the 
whole particle distribution within the observational cylinder, before and after 
the inclusions of the observational biases, respectively. Since the two simula- 
tions are started with the same initial random numbers and the two clusters 
lie almost at the same position, they are originated by the same waves and any 
difference in their morphology is only due to the difference in the correspond- 
ing cosmologies. The cluster in ACDM0.3 has a much better defined shape than 
that in TCDM. The latter, even after the interloper removal, does not look 
in projection like a well defined structure. As a consequence, the membership 
assignment is not very efficient in identifying genuine cluster members and the 
resulting velocity dispersion profile does not recover the true a^ so efficiently 
hke for the ACDM0.3 cluster. 

The results of this analysis can be summarized into two main points. 

(a) Although the scatter between the intrinsic and the observed a^ is in gen- 
eral rather large, the difference between the corresponding N{> o"„) is not 
dramatic, especially for o"„^ 800kms~^. This is just the consequence of the 
roughly symmetric distribution of overestimates and underestimates of a^. 
Discrepancies in the high a„ tail are due to few clusters whose a^ ought to 
be overestimated by 200-300 km s~^, because of the presence of unremoved 
interlopers. A typical example of this occurrence is just provided by the 
third cluster in Figure 2. 

(b) Differences between intrinsic and observed A^(> a„) may however not 
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Fig. 5. The same as in Figure 3, but for CHDM, OCDMo,4 and OCDMo.e simula- 
tions. 

be negligible. Furthermore, no a priori recipe exists, which could allow to 
recover the correct CVDF from the observed one, the difference being non- 
trivially model dependent. This casts doubts on the reliability of detailed 
comparisons between the abundance of clusters, as inferred from their inter- 
nal velocity dispersion, and predictions of DM models based on analytical 
approaches, like the Press & Schechter (1974) one, which can not include 
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any observational effect. 
4-3 Comparing with DM models 



Based on the results obtained in the previous sections, we compare in Figure 7 
the CVDF for the observational sample and for the simulations. In each panel, 
the shaded band represents the result of our analysis of the F96 sample (cf. 
Figure 1 and Section 4.1). For each model, results are obtained by averaging 
over the three available realizations (exept for CHDM). 

As for the SCDM model, as expected it largely overproduces clusters at any 
cr„. In order to check for the effect of a possible velocity bias on this model, we 
resorted to the relation by = 0.7(r/0.5 /i~^Mpc)°-^ (6^,: ratio between the ve- 
locity dispersions of "galaxies" and dark matter particles) obtained by Evrard 
et al. (1994) from hydrodynamical cluster simulations of a CDM model with 
as = 0.6. Since these authors found the above relation to be almost indepen- 
dent of the evolutionary stage, we adopt it also for our larger as output (cf. 
Table 1). The effect of such a velocity bias is shown with the dotted curve in 
the upper left panel; even in this case the resulting A^(> a^) is much larger 
than the observational one. 

From the results reported in lower left panel, it turns out that lowering the 
primordial spectral index to n = 0.85 is not enough to bring a r2o = 1 CDM 
model into agreement with observations, unless one introduces a substantial 
amount of gravitational wave contribution to the CMB anisotropy. Even tak- 
ing T/S = 7(1 — n) from power-law inflation {as = 0.67), the CVDF remains 
significantly larger than that of the observational sample. Only further low- 
ering the normalization to as = 0.51 leads into agreement with real data, 
thus in accordance with the X-ray based results (see, e.g., E96; Pen 1996). 
In their analysis of the large-scale velocity fields, Zaroubi et al. (1997) found 
that TCDM with n ~ 0.85 is the best fit for ^o = 1, ^i = 0.5 CDM models 
only if T/S = (see, however, Borgani et al. 1997b). A possible way out may 
be allowing for a substantially larger baryon fraction (see, e.g.. White et al. 
1996). This would have the effect of suppressing fluctuations on the cluster 
mass scale, while leaving the spectrum almost unchanged at the larger scales 
probed by velocity flelds. 

As for the CHDM model, it turns out to marginally overproduce clusters. A 
better agreements can be achieved either by increasing the hot component 
or by tilting the primordial spectral index. In both cases this model, which 
is already marginally consistent with observations of high-redshift {z = 3-4) 
damped Lya systems, would have even harder time with the galaxy forma- 
tion timing (see, e.g., Borgani et al. 1996, and references therein). A way to 
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Fig. 6. Comparison of the role of projection effects on a cluster identified in ACDMq.s 
(upper panels) and TCDM with ag = 0.51 (lower panels) simulations. Note that 
the two clusters have almost the same positions in the two simulations based on 
the same set of initial random numbers. Left panels show in projection the whole 
distribution of simulation particles within the observational cylinder. Central and 
right panels are the same as for plots in the first and fourth lines of Fig. 2. 

overcome this problem would be to share the hot component between more 
than one z/ species; in this case, the increase in the neutrino free streaming 
suppress fluctuations on the cluster mass scale, without significantly affecting 
the fluctuation power on the galaxy mass scale. A model with Q^ = 0.2 and 
two massive neutrinos (Primack et al. 1995; Primack 1996) has been found 
to satisy at the same time the constraints from the cluster abundance and 
high-redshift objects. 

In any case, it is interesting to note that the CHDM model, for which as = 
0.77, has a CVDF which is comparable to (or even smaller than) that of 
TCDM with as = 0.67, even though fio = 1 for both models. This is the 
consequence of the presence of the neutrino component, which acts like a sort 
of background superimposed on the more clustered cold component. Its effect 
is that of slowing down the dynamical clustering evolution, so that more time 
is required for a CHDM model to develop the same small-scale velocities of a 
purely cold model. This means that a larger as value is required for a CHDM 
model to produce the same cluster abundance as an fio = 1 CDM model. 

As for the ACDM models, the Qq = 0.3 case is quite consistent with the 
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Fig. 7. Comparisons between CVDF's for real data (dashed band) and simulations. 
Simulation results correspond to the average over the three random realization of 
each model (except for CHDM, for which only one realization is available). 

observational A^(> cr^), thus confirming the result based on X-ray data about 
the capability of this model to produce the correct cluster abundance (see, 
e.g., E96). The ACDMq.s model, instead, turns out to overproduce clusters, as 
expected on the ground of its large Ug value. As for OCDM, the Qq = 0.4 model 
is found to underproduce clusters, while increasing the density parameter to 
Qq = 0.6 turns into a marginal cluster overproduction. 



5 Conclusions 



In this paper we addressed the question concerning the reliability of the dis- 
tribution of cluster internal velocity dispersions as a diagnostic for cosmo- 
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logical models. In order to properly address this issue we decided: (a) to 
run simulations for a variety of models, so as to verity the discriminative 
power of A^(> a^); (b) to analyse simulated clusters in the same way as a 
reference observational sample, originally considered by Fadda et al. (1996; 
F96), consisting of 153 clusters and complete in velocity dispersion down to 
(^v,iim - 650kms~^. 

Our main results can be summarized as follows. 

(a) Projection effects and the limited number of galaxy redshifts per cluster 
can heavily pollute the recovering of the correct cr^. The cluster member 
selection procedure leads to the wrong identification of several interlopers as 
genuine members. This effect is more pronounced for low-a^, {^ 600kms~^) 
objects. On the other hand, high-o"^ clusters, which look in projection like 
better defined structures, display an improved recovering of the true a^. 

(b) The resulting effect is in general that of overestimating the CVDF, espe- 
cially in the high-cr^, tail. This is due to few small, Iow-ct^, clusters, whose 
velocity dispersion is heavily boosted by projection effects. Furthermore, the 
amount of the overestimate is non-trivially model dependent: simulations in 
which particles last more within virialized structures generate clusters with 
a sharper profile, even after projection, so as to make interloper removal an 
easier task. 

(c) As for the comparison with data, our results substantially agree with the 
(jg-fio scaling for CDM models, based on the Press-Schechter approach to 
the X-ray temperature function (e.g., Viana & Liddle 1996; Eke et al. 1996; 
Pen 1996). We confirm that a fig = 1 CDM model requires as — 0.5. The 
ACDM model with Qq = 0.3 produces the correct N{> ay) with as = 1.1. 
This value is somewhat larger than that predicted by Eke et al. (1996; cf. 
Table 1), but rather consistent with the scaling provided by Pen (1996). 
As for the CHDM model, it overproduces clusters, although not by a large 
amount. In any case, the presence of the neutrino background has the effects 
of slowing down the development of non-linear structures. As a consequence, 
CHDM models are allowed to have a larger as value to provide the same 
cluster abundance as an f^o = 1 CDM model (see also Jing & Fang 1994 
and Walter & Klypin 1996). This is the reason why the CVDF for CHDM 
with as = 0.77 is smaller than that of TCDM with erg = 0.67. 

Such constraints on ag from cluster abundance can be compared with 
those from large-scale velocity fields. For instance, Zaroubi et al. (1997) 
found for CDM models that ag^lj-^ = 0.88 ± 0.15 (90% c.l.) is predicted 
by a maximum likelyhood approach to the peculiar velocities of the Mark 
111 sample. Instead, Borgani et al. (1997b) analysed the sample of cluster 
pecuhar velocities by Giovanelli et al. (1997) and derived asfi^'^ = 0.5 ±0.2 
(90% C.L), thus in better agreement with cluser abundance. 

Based on such results, we conclude that cluster internal velocity dispersions 
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represent a stringent cosmological constraint, provided that observational bi- 
ases, like projection and sampling effects, are carefully accounted for (cf. also 
van Haarlem et al. 1997). It is also clear that further advancements both from 
the theoretical and the observational sides would be required in order to fully 
exploit the potential of cluster velocity dispersions. For instance, the availabil- 
ity of high-resolution simulations, also including gas dynamics, should settle 
the question of the velocity bias, whose amount should be rather model de- 
pendent. 

Furthermore, improved observational strategies should eliminate some of the 
ambiguities, which are still present in the data sets. The possibility of measur- 
ing many more galaxy redshifts per cluster will allow a better understanding 
of the internal cluster dynamics and of the role of substructures. Extensions 
of the today available samples in the following two directions would provide 
even more stringent constraints on DM models. From the one hand, realis- 
ing a shallower volume-limited {z^ 0.1) survey, including at least down to 
R = clusters, would extend to smaller a^ the range where the CVDF is 
accurately sampled. From the other hand, realising a deeper {z^ 0.5) survey 
of rich clusters would allow to discriminate between those models that, even 
though providing the correct A^(> a^) at z ~ 0, have different evolutionaly 
patterns (e.g., Jing & Fang 1993; Crone & Geller 1994). Attempts in this di- 
rection have been already pursued with the CNOCC cluster survey (see, e.g., 
Carlberg et al. 1996, and references therein), but only for a limited number of 
clusters. In any case, it is clear that as deeper cluster surveys are considered, a 
careful treatment of projection biases becomes a more and more delicate issue. 
Multi-object spectrographs of the new generation will be needed to measure 
~ 10^ galaxy redshifts for each cluster in one shot, thus rendering feasible in 
the near future a substantial enlargment of the today available data sets. 
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